Transport properties of the quark-gluon plasma from lattice QCD 



Harvey B. Meyer 

Center for Theoretical Physics 
Massachusetts Institute of Technology 
77 Massachusetts Ave 
Cambridge, MA 02139, U.S.A. 



Abstract 

I review the progress made in extracting transport properties of the quark-gluon plasma from lat- 
tice QCD simulations. The information on shear and bulk viscosity, the "low-energy constants" 
of hydrodynamics, is encoded in the retarded correlators of T^, the energy-momentum tensor. 
Euclidean correlators, computable on the lattice, are related to the retarded correlators by an in- 
tegral transform. The most promising strategy to extract shear and bulk viscosity is to study the 
shear and sound channel correlators where the hydrodynamic modes dominate. I present prelim- 
inary results from a comprehensive study of the gluonic plasma between 0.95 T c and 4.07V 



1. Introduction 

The phenomenology of heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC) 
has revealed unexpected properties of the quark gluon plasma. Hydrodynamics calculations U 
successfully described the distribution of produced particles in these collisions [2]. This early 
agreement between ideal hydrodynamics and experiment has been refined in recent times. The 
dissipative effects of shear viscosity 77 have been included in (2+l)D hydrodynamics calcula- 
tions and the sensitivity to initial conditions quantitatively estimated [4] for the first time. 
Experimentally, the elliptic flow observable V2, which is sensitive to the value of 77 in units of 
entropy density s, is now corrected for non-medium-generated two-particle correlations 0. The 
conclusion that 77/ s must be much smaller than unity has so far withstood these refinements of 
heavy-ion phenomenology |6 1. On the theory side, it is therefore important to compute the QCD 
shear viscosity from first principles in order to complete the picture. Furthermore, since the 
heavy ion collision program at LHC will probe the quark-gluon plasma at temperatures about a 
factor two higher Q, it is crucial to predict the shear viscosity at ~ 3T C , and to relate it to the 
size of elliptic flow, before experimental data is available. 

In this talk I present lattice calculations of the thermal correlators of the energy-momentum 
tensor (EMT) in the Euclidean SU(3) pure gauge theory (i.e. quarkless QCD), and discuss meth- 
ods to extract the shear and bulk viscosity from them. Computationally, the calculation is chal- 
lenging enough without the inclusion of dynamical quarks, and physically, the static properties 
of the QGP, normalized by the number of degrees of freedom, do not depend sensitively on the 
flavor content away from T c |8|. In perturbation theory [9|, the ratio of shear viscosity to entropy 
density is O(^) and there is only about 30% difference between the pure gauge theory and full 
QCD 1 10] at a fixed value of a s (tj/s is smaller in the pure gauge theory). The bulk viscosity £ is 
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O(ffj) and predicted to be much smaller than the shear viscosity in perturbation theory. As a rule 
of thumb, T]/s « 1.0 g) and £/t] ~ 10~ 3 Q3] at a s = 0.25. 

Since the pioneering calculation [12], there have been only few attempts to calculate the 
viscosities of the pure gauge theory on the lattice. Nakamura and Sakai |[T3l performed the first 
calculations at N T = (aT)~ l = 8, where a denotes the lattice spacing. About two years ago, the 
accuracy of the Euclidean correlators improved significantly [ 14, 15] thanks to high statistics and 
a more efficient two-level algorithm |[T6l [T71 . Here I will give an update on the progress made 
in constraining dynamical properties of the QGP such as the viscosities. I will not describe the 
calculations of the electric conductivity [18| or heavy quark diffusion constant lfT9ll . except for 
saying that they share many features with the viscosity studies. 



2. Hydrodynamics and energy-momentum tensor correlators 

From the modern point of view, hydrodynamics is an effective theory that describes the slow, 
long-wavelength motion of a fluid. The central object is the energy momentum tensor, a symmet- 
ric rank-two tensor, whose conservation is expressed by the four continuity equations d^T^' = 0. 
The component Tqq is the energy density, Tot is the momentum density, which coincides with 
the energy flux, and the spatial components are the momentum fluxes. From the micro- 
scopic point of view, the matrix elements of T^ y , viewed as a quantum operator, between any two 
on-shell states satisfy d M (V\T^(x)\<S>) = 0. 

The two macroscopic modes of fluid motion near equilibrium are the shear and sound modes. 
The former corresponds to the transverse diffusion of momentum. Indeed a small perturbation 
7*03 = To3(t,x{) of the fluid around equilibrium gives rise to shear flow and satisfies the diffusion 
equation 

d, Tw(t, xi) - Dd\ T 03 (t, xi) = , (shear mode) (1) 

where the diffusion coefficient D is proportional to the shear viscosity r\, D = ^ (e is the 
energy density and p the pressure). Secondly, a sound wave propagating in the z-direction with 
wavelength A = 2n/k is damped according to 

T 03 (t, k) oc e -\(h+ok 2 tKe+ P ) (sQund mode) _ (2) 

Both the shear viscosity rj and bulk viscosity f thus contribute to the damping of sound waves. 

The quantities computed on the lattice are correlators of the energy-momentum tensor that 
depend on Euclidean time xo and spatial momentum q, 

C^pAxo, q) = L 5 J dH e iqx {T^xq, x)T pa (0)) , L =l/T. (3) 

The so-called spectral function piu, q) is defined as (— n x) the imaginary part of the Fourier- 
space, Minkowski-time retarded correlator Gfi(<jj,q). I leave the temperature dependence of 
the functions C and p implicit. For fi — p and v = <x, it obeys a positivity condition, in our 
conventions i"p(aj,p)/a> > where n is the number of time components among the four indices, 
and is odd in to, p(-u>, q) = ~p(o), q). Via a Kubo-Martin Schwinger relation, this spectral 
function is related to the corresponding Euclidean correlator (see [ 19 1), 

r°° cosh oj{\L - x ) 

C^pa-(xo,q) = L Pp V ,pJco,q) — — j da>. (4) 

Jo sinh jO)Lo 
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A wealth of information is encoded in the spectral functions. In particular, the shear and bulk 
viscosities are given by Kubo formulas (see [20 1 for a derivation), 

,™ r ,. Pl3,13(w,q) 4 P33,33(^,q) 

77(7) = 7r lim hm , Z£(T ) + 7](T ) = n limhm . (5) 

«->0q->0 CO ' «->0q-<0 CO 

They also determine the static structure of the plasma [21 , 22], 

CWr(r, D = (T MV (x , r) r ptr (* , 0)) = lim f £\ f dco e^ 9 ^^^ . (6) 

e->0j (27r) j Jo tanhco/zT 

In the following, I summarize the known analytic properties of the spectral functions. The goal 
will be to exploit these properties to help us solve the integral equation |4]) for the spectral func- 
tion, given the numerically calculated Euclidean correlators. 

The conservation of the EMT, c^7), v = 0, implies in particular, for q = qe^, 

w 4 poo,oo(w,q) = -0) 2 q 2 p 03fi3 (a),q) = <7 4 P33,33(«,qX -u? p om {o>, q) = g 2 pi 3 , I3 (w,q) . (7) 

The first set of spectral functions corresponds to the sound channel and the second to the shear 
channel [20]. I therefore use the notation p sn d(w, q) = -po3,03(w,q) and p s h(w, q) = -poi,oi(w,q). 

Kovtun and Starinets Il23l analyzed the general tensor structure of EMT correlators. I briefly 
summarize their results. For a generic relativistic quantum field theory in four dimensions, there 
are five independent functions of (co, q) that determine all thermal correlators of the EMT. At 
T = 0, this number is reduced to twcQ[24]. A physically motivated choice of five independent 
functions is: two functions corresponding to the two hydrodynamics modes, the shear mode (Gi 
in the notation of [23 1) and sound mode (2G2 + 3C/J; one function associated with the {T^Tn) 
correlator and momentum in the z direction (G3); finally, the correlators (T^Tyy) and (T^Too) 
(respectively proportional to 2C7- + Cl and Cj + 2Cl) determine the last two functions. As an 
application, it is worth noting that if 9 q is the polar angle of q, 



d(D {p S nd(w, q) cos 2 q + p sh (w, q) sin 2 q } 



cosh<y(jLo - *o) 
sinh ^o>Lq 



(8) 



In the future this will allow us to exploit also momenta that are not aligned along a lattice axis. 

At T — 0, G\ — G2 — G3 and Cj = Cl, and they become functions of {u> 2 - q 2 ). As a 
corollary, the extent to which these relations are obeyed at T > provides us with a way to 
evaluate the importance of thermal effects on the correlation functions. 

As a low-energy theory, hydrodynamics predicts the small k = (o>, q) behavior of spectral 
functions in terms of a few 'low-energy constants'. A first-order hydrodynamic expression for 
the shear and sound spectral functions was derived in [20]. For q = qej,, 

Psh(^, q) u>, q ->o rj_ <f_ 

co n a) 2 + (i]q 2 /(e + P)) 2 ' 

Psnd(w, q) oj,q^0 r s (e + P) q 2 LJ 2 



n (to 2 - v)q 2 ) 2 + (T s coq 2 ) 2 



(10) 



'For a conformal field theory, the number of independent correlators is three at T > 0, and just one at T = 0. 
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Figure 1: Left: Euclidean sound channel correlator. Scaling of CPU-cost with N r in the deconfined phase. 



I introduced the speed of sound v s and the sound attenuation length F s = (|/7 + £)/(e + Baier 
et al. [25 1 derived the dispersion relation of the sound pole to next-to-leading order accuracy for 
a conformal theory: 



s (q)q, 



s (q) 



l + 



1 \Tn 



4v? 



0(q 4 ) 



(11) 



Here Tn is the relaxation time for shear stress. 

The Euclidean correlators and spectral functions were calculated to leading order at weak 
coupling in 1 26 1 . The correlators are also known to next-to-leading order at zero temperature [27, 
[24). Expressing the latter results in terms of a shear channel spectral function, 



P13,13(W,0) 



d A a> 
10(4^) 2 



1 



5a s N c 
9n 



(12) 



Son and Romatschke [28| derived two sum rules for the spectral functions. The 'bulk' sum 
rule can be reexpressed in terms of sound and shear channel spectral functions, 



f 

Jo 



dot 



OJ 



P33,33(W,0) -pi3,i 3 (w, 0) 



T-0 



From a practical point of view, this provides an additional constraint on the spectral functions, 
albeit at the price of having to determine simultaneously the T-0 spectral functions. Because 
the symmetry group is larger at T = 0, the vacuum spectral functions can however be more 
strongly constrained by lattice data (in addition to having a longer Euclidean time extent). 

In a class of strongly coupled gauge theories at large-rV c , exact results have been obtained 
for the transport coefficients. The most important result is that the shear viscosity assumes a 
universal value in a large class of theories, rj/s = l/4n [29|. The spectral functions for the EMT 
have also been obtained by numerical integration in the case of the N — 4 SYM theory |30|. 
The result is compared to the hydrodynamic functional form (lOi in (|21|, Fig. 3). Such a 
comparison teaches us up to what frequency and momentum the hydrodynamic spectral function 
remains a good approximation to the exact result. It turns out that up to w,^ w jtT, the sound 
channel spectral function is well described by hydrodynamics. This sort of statement necessarily 
depends on the details of the microscopic theory; and in a weakly coupled theory, ones expects 
the validity of hydrodynamics to be limited to much greater wavelengths. Nevertheless I will use 
this observation as a guideline for the gluon plasma at T < 4T C . A similar logic has also been 
applied in studies of Mach cones in the plasma BTl . 
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Figure 2: Sound channel spectral function at 2.3T C . In the Table, lattice results are compared to free gluons [33] and to 
the strongly coupled SYM results ['34 25]. Stat, errors only are given. We expect ( to be negligible at these temperatures. 



3. Lattice calculation 

I use Monte-Carlo simulations of the anisotropic Wilson action for SU(3) gauge theory with 
a fixed anisotropy of £ = 2, meaning that the temporal lattice spacing a T is half the length of 
the spatial lattice spacing a, T . See l32l for details of the action. The lattice spacing is related to 
the bare coupling (an input parameter of the simulation appearing in the lattice action) through 

~ 1/ log(l/a -A). Due to the loss of continuous translation invariance on the lattice, the dis- 
cretized EMT is not protected against a non-trivial go-dependent renormalization. Fully calibrat- 
ing the EMT on the anisotropic lattice remains a numerical challenge. While To* is normalized 
multiplicatively and 7oo and can be calibrated by using thermodynamic expectation values, 
the spatial components present the biggest challenge. 

A typical example of a Euclidean correlator computed on the lattice is shown in Fig. [T] It 
displays the correlator of the energy density Too, at various non-vanishing spatial momenta. The 
lattice size is 16 x 48 3 and the temperature about 2.37V The typical size of the error bars is 
3%. Obtaining the correlators to that accuracy is computationally demanding. Figure [T] also 
shows how the cost of the computation scales with N T . Different temperatures lie approximately 
on a smooth curve. The dashed line corresponds to the expectation for a conformal theory, 
tcpu 00 N^N^.. For a fixed aspect ratio N T /No-, the cost therefore scales as N^. This is in contrast 
with the expectation of an N^N^. scaling with a one-level algorithm. 

The discretization errors affecting the Euclidean correlators calculated at finite lattice spacing 
were studied in perturbation theory in [35 1. The study motivated in particular the choice of 
anisotropy flo-/a T = f = 2. Further it was shown that it is advantageous to remove the treelevel 
cutoff effects from the lattice correlators. I apply this technique to all the data shown in these 
proceedings. Finally, at treelevel the q = correlator of Tot is much flatter as a function of Xq 
than the correlator of Too (they become exactly independent of Xq m me continuum limit). This, 
and the fact that Tm is normalized multiplicatively even on the anisotropic lattice, means that the 
correlators of the momentum density are the most accurately determined of all. 

Relations |7]i imply that in the sound channel, the numerically determined correlators Cqo.oo^ 
Coo,03> ^03,03, C33_o3 and C3333 are all related to the same spectral function p sn( j modulo factors 
of (u>/q). Therefore a global analysis of these five correlators is the best way to constrain p sn d, 
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Figure 3: Left: Coo.oo(*o, l) above and near T c , where evidence for the damping of sound waves is seen. Right: test of 
a Ward identity at finite lattice spacing. 



provided the Euclidean correlators are consistent with the Ward identities. This has to be the case 
in the continuum limit, but presently analyses are still based on data at finite lattice spacing. It is 
therefore essential to check in at least some cases that the Ward identities are satisfied to an ac- 
curacy comparable to the statistical errors. Conversely, these identities can also be used to deter- 
mine some of the normalization factors of T MV . Figure ^ compares the correlators Co3,o3(*o>q) 
and Coo,33(*o> q) f° r several values of Xq and q = qej, (they are exactly equal in the contin- 
uum). While the normalization of Tib is easily determined by requiring Co3 ; o3(*o,0) = s/r 3 , 
the operator ^3 requires in total four unknown normalization factors. Two of them are fixed by 
thermodynamics, and two have been fitted so that the two correlators on Fig. Q agree. Since 
there are O(20) points where the correlators are compared, this represents a non-trivial check. 

In the following I present an analysis of the sound channel correlators on a 16 x 48 3 lattice 
for xq > 1 /AT (corresponding to xo/a T > 4) and for q < jiT . I assume for now, motivated by the 
inspection of the strongly coupled N - A SYM spectral functions, that hydrodynamics describes 
the low-frequency part of the spectral function up to that momentum. 

I adopt a fit ansatz motivated by perturbation theory at high frequencies (phighX and hydrody- 
namics at low energies (pi ow ). Specifically, p snd = p low + p med + p high , where 



pi ow ((y, q, T) 2F S (e + P) 10 q 1 + <xi u> 



(14) 



tanh(w/2J) n ( w 2 _ v ^ q ) q 2f + (T^ 2 ) 2 1 + o"2« 2 

p me d(a>,q,T) 2 2 2 0) £cr 

= a) q tanh ( — )—z ~ z , (15) 

tanh(w/2r) H 2T V 2 + (w 2 - q 2 - M 2 ) 2 V ; 

Phigh(«, <7, T) 22 ,2 w 2d A 

: tz j^fz = o) q tanh — — r 2 < ^ 

tanh(w/2r) 2T 15(47r) 2 



where v s (q ) is given by Eq. 1 1 Some remarks are in order. Perturbation theory and the operator 
product expansion can be used to systematically improve the knowledge of phigh- Secondly, a full 
second-order hydrodynamics parametrization of the spectral function at low-frequencies would 
significantly improve our understanding of the accuracy on the transport coefficients that can be 
expected from such as global fit. Thirdly, the region a> = 0(T) is the one where there is least 
theoretical guidance. Given the present amount and quality of the data, the simple ansatz (JT3J 
used here is a reasonable choice (other functional forms could be used to test the sensitivity of 
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the transport coefficients to this particular choice). As the quality of the data improves, the goal 
is to treat this part of the spectral function (p me d) in a more systematic way, either by expanding 
it in a basis of orthogonal functions [ 14] or by using the Maximum Entropy Method [36 1. 

In the present specific example I fit 7 parameters to a total of 48 data points. The parameters 
are F s , <j\ , cr 2 , Tn, C, cr, M. Figure ^ displays the reconstructed sound spectral function based on 
lattice data at 2.3T C , and the results for the fit parameters of interest are given there too. 

There has been significant interest in the behavior of bulk viscosity near the phase transi- 
tion ll37l l38l 1391 . Figure ^ left panel) displays the Coo.oo correlator as a function of q 2 . Re- 
calling the relation between the Euclidean correlator and the spectral function (Eq. [4]) as well 
as the positivity property of the latter, this Figure shows that for T = 1.0 17^ and q = 7tT/4, 

~cPr* Jo^^ 77P»nd(fc>> Q) c v /cIaT 3 ~ 13.6(1.9) (c v is the specific heat). The sound peak has 
therefore quasi disappeared even for a sound wavelength of A — 8/T w 5.8fm. I conclude that 
for all practical purposes, the gluonic medium does not support sound waves near T c . 



4. Conclusion 

I have described the calculation of transport properties of the quark gluon plasma on the lat- 
tice. For the time being the calculations are performed in the purely gluonic plasma in order to 
reduce the computational cost. The primary quantities computed are the Euclidean-time depen- 
dent correlators of the energy-momentum tensor. Based on these the spectral functions can be 
constrained, particularly with the help of information on their functional form. Hydrodynamics 
and perturbation theory describe (in opposite regimes) both the spatial momentum and frequency 
dependence of the correlators. Correlators with finite spatial momentum are then useful, because 
within either the shear or the sound channel they are interrelated by the energy and momentum 
conservation equations. Different correlators can thus be used to over-constrain the shear and 
sound spectral functions. Altogether, it is now possible to include 50-100 data points in a global 
analysis of the sound channel, and, as illustrated by Eq. [8] even more points per channel can be 
used if the shear and sound channels are analyzed simultaneously. 

My current best guess for the viscosity of the hot matter that will be created at LHC uses 
the lattice result rj/s as 0.26 for the purely gluonic plasma at 2.3T C , and multiplies it by the 
perturbative (AMY) ratio of rj/s in full QCD and in the pure gauge theory |9|, 



[ll/s]QCP ~ [T]/s]GP,lattice 



QGP 



* 0.40. (17) 

AMY 



I have presented direct evidence for the disappearance of sound waves as collective excita- 
tions of the plasma in the vicinity of the deconfining phase transition. This effect leaves a strong 
signature on the Euclidean correlators in the pure gauge theory, and it would be very interesting 
to know how strong it remains in full QCD, where the transition to the QGP is a crossover. As the 
chiral critical point is approached at finite baryon density (assuming it exists), the suppression 
of sound waves becomes stronger. This has recently been studied in detail by a combination of 
hydrodynamics and the theory of critical phenomena [40 1, and even been proposed as a possible 
signature of the chiral critical point. 

I thank the organizers of the Quark Matter 2009 conference for a very enjoyable and stimulat- 
ing physics event and for the opportunity to present this work. Lattice computations were carried 
out on facilities of the USQCD Collaboration, which are funded by the Office of Science of the 
U.S. Department of Energy, as well as on the Blue Gene L rack and the desktop machines of the 
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Laboratory for Nuclear Science at M.I.T. This work was supported in part by funds provided by 
the U.S. Department of Energy under cooperative research agreement DE-FG02-94ER40818. 
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